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Abstract 



Oh! 

r~| ! Self-consistent approximations in terms of fully dressed propagators provide 

. ,_, a simple expression for the entropy of an ultrarelativistic plasma, which iso- 

H lates the contribution of the elementary excitations as a leading contribution. 

Further approximations, whose validity is checked on a soluble model involv- 
ing a scalar field, allow us to calculate the entropy of the QCD plasma. We 
obtain an accurate description of lattice data for purely gluonic QCD, down 

to temperatures of about twice the transition temperature. 
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The properties of the high temperature phase of QCD cannot be easily calculated using 
perturbation theory, in spite of the fact that the gauge coupling g is small if the temperature 
T is sufficiently high. This is evidenced in particular by the poor convergence properties of 
the perturbative series |I[. 

Lattice results, which show that the ideal gas limit is approached as T becomes large, 
can be accounted for reasonably well by phenomenological fits involving massive quasi- 
particles |@,|1. Although the quasiparticle picture suggested by such fits is a rather crude 
representation of the actual physics of non-abelian gauge theories, it supports the idea that 
one should be able to give an accurate description of the thermodynamics of the QCD plasma 
in terms of its elementary excitations. 

It is worth emphasizing at this stage that, among the relevant degrees of freedom, the 
soft collective ones, with momenta of order gT, are clearly non-perturbative. Although their 
leading order contribution oc g^ to the pressure can be easily isolated 0], it does not make 
much physical sense to regard this contribution as a genuine perturbative correction. Recent 
investigations |§|-0| indicate indeed that trying to represent this contribution by a truncated 
polynomial in g is not appropriate. 

In order to carry out a more complete calculation, we shall use techniques which allow 
systematic rearrangements of the perturbative expansion, avoiding double countings. We 
shall rely in particular on self-consistent approximations which provide a simple expression 
for the entropy, isolating the contribution of the elementary excitations as a leading contri- 
bution. We show that this entropy formula can be used to get a good estimate of the QCD 
entropy at high temperature. The results that we have obtained so far are quite encouraging 
and give us the hope that an analytical control of the high temperature phase of QCD is 
within reach. 

We shall first discuss scalar field theories with gip^ + g'^i^^ interactions, with the double 
purpose of presenting the general framework and of checking approximations which will be 
used later for QCD. 

The thermodynamic potential Vt = —PV of the scalar field can be written as the following 



functional of the full propagator D PQ : 

pn[D] = -logZ = iTr logD-i-iTrHD + $[£)] (1) 

where Tr denotes the trace in configuration space, (3 = 1/T, and $[-D] is the sum of the 
2-particle-irreducible "skeleton" diagrams 

-$[D]= 1/12^— ^ +1/8 r'^^ +1/48^3 +... (2) 

The self energy 11 = D^^ — Dq^, where Dq is the bare propagator, is related to $[-D] by 

6(!?[D]/6D = in. (3) 

An important property of the functional f2[-D], which is easily verified using @, is that it is 
stationary under variations of D: 

6n[D]/6D = 0. (4) 

Self-consistent ( "^-derivable" ) ^ approximations, i.e., approximations which preserve this 
property, are obtained by selecting a class of skeletons in $[-D] and calculating 11 from Eq. (^ 
above. 

The stationarity of Q[D] has an interesting consequence for the entropy. Because of 
Eq. (H) the temperature derivative of the spectral density in the dressed propagator cancels 
out in the entropy density 

S = -d{Q/V)/dT (5) 



and one obtains 10.11 



d'^k dn{uj 



J 27iY dT ^ ^ ' ^ 



(27r)4 dT 

f d^k dn{u. 

(27r)4 dT 
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with 



9p W«^R,ntofl = (7) 

dT D J (27r)4 dT ^ ' 
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up to terms that are of loop-order 3 or higher. Thus, in contrast to Q, where $ contributes 
already to order g"^ in perturbation theory, Eq. (5) with 5' = is perturbatively correct 
up to, and including, order g^ [|12|. The first two terms in Eq. (^) represent essentially the 
entropy of "independent quasiparticles" , while S' may be viewed as the residual interactions 
among these quasiparticles |TT[ . 



Besides this important simplification, S, in contrast to the pressure, has the advantage 
of manifest ultra-violet finiteness, since dn/dT vanishes exponentially for both u -^ ±oo. 
Moreover, any multiplicative renormalization D -^ ZD, U -^ Z~^Il with real Z drops out 
from Eq. (|). 

We now focus on the self-consistent approximation obtained in g'^(f^ theory where only 
the second diagram in Eq. (H) for $ is kept. Then Imll = and Rell = m^ = const., and 
Eq. (Bl) reduces to 



5=-/igj^Imlog(*^-.^ + ™=) 
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A fully self-consistent determination of m corresponds to solving the gap equation 

m' = 12^2 jAA_n(u;)e{u)S{u^ - ^ - m"). (9) 

When the solution of this equation is inserted in Eq. (H), the entropy obtained coincides 
with that of the exact solution of a scalar 0(A^)-model in the limit of A^ ^ oo |[T3|,|7[|. Note 



that in contrast to Eq. (||), the gap equation is ultraviolet divergent (for uj -^ — oo), and 
requires renormalization 0, affecting m? at perturbative order g'^ and beyond. 

In view of the subsequent application to QCD, where a fully self-consistent determina- 
tion of the gluonic self-energy seems prohibitively difficult, we consider now perturbative 
approximations. Our goal is to obtain approximate expressions for the self energy which 
allow us to reproduce the perturbative result for the entropy when expanded to order g^. 
We emphasize that our final results for the entropy are non-perturbative, and not limited 
to a truncated polynomial in g. What we are testing here is the quality of approximations 



which preserve self-consistency up to order g^ at least, which is what we shall be able to do 
in QCD. 

As first approximation we consider the leading contribution to the self-energy at high 
temperature, the so-called hard thermal loop (HTL) [|1^], and refer to Eq. (^ with this 
restriction as iShtl- For the (p^ theory, we have simply IT ^ 11 = -m? = g'^T'^. When inserted 
into Eq. (|^), this yields the correct result for the leading-order interaction term g^T^ in the 
entropy. 

On the other hand, the order g^ contribution contained in iShtl turns out to be too 
small by a factor of 4 when compared to the well- know perturbative result 0]. This is 
corrected by including the next-to-leading order (NLO) term in the thermal mass through 
resummed perturbation theory, m^ = m^ + 6m'^ = g'^T'^ — -g^T^ 0]. So, when compared to 
conventional resummed perturbation theory, the order-^f^ term arises in an unusual manner: 
whereas in the former the entire plasmon effect comes from the infrared regime, in Eq. (^ 
an even larger contribution comes indirectly from the infrared through corrections to the 
dispersion laws relevant at hard momenta (the T^ term in Eq. (^ comes entirely from hard 
A; ~ T). This may be understood as a consequence of the requirement of self-consistency: 
recall that Eq. (^ relies on the stationarity of the thermodynamic potential, and this has 
to be maintained at the order of interest. 

At large coupling the NLO result for m^ inevitably turns negative. This can be avoided 
by taking instead the perturbatively equivalent form rfi^ + 6m'^ = g'^T'^ /[I + 3(7/7r], which is 
monotonous in g, giving a very good approximation to the solutions of Eq. (H) up to (7 > 1. 

In Fig. 1 we compare the various approximations numerically with the exact entropy 
(full line), normalized to their ideal-gas values (SB), as functions of the renormalized cou- 
pling in the MS-scheme. In contrast to the full result, the perturbative approximations are 
renormalization-scale dependent. As in Ref. ^ we consider the effect of varying the scale /i, 
here in the range /2 = (| . . . 2) x 2-nT. The lower and upper dark-gray bands correspond to 
conventional perturbative results for S/ S^b up to order g"^ and g^, resp., the medium-gray 
bands to the result of HTL- and NLO-resummations of the two-loop entropy. The latter. 



which clearly represent a substantial improvement over the former, are the approximations 
that we shall now implement for QCD. 
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FIG. 1. Comparison of perturbative and HTL-improved approximations to the entropy in the 

large- A^ scalar 0(A^)-model. See text for detailed explanations. 

The analog of Eq. @ in purely gluonic QCD and in Coulomb gauge reads 
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with A^„ 



Ar2 



1 = 8 for SU(3). The (spatially) transverse (Ht) and longitudinal (Hl) 



structure functions will be specified below. Here we have assumed that the gluon self-energy 
is transverse with respect to the four-momentum, and that there are no contributions from 
the ghosts, which turns out to be justified in the approximations that we shall be interested 
in p|. 

The order g"^ contribution to the entropy is easily extracted from Eq. (p!oD: 

d'^k On 
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Here we have used the fact that the integral is dominated by hard momenta and that the 
transverse quasiparticles have the asymptotic thermal mass: 



mi, = IIt{uj' = e) = g'NT'/6. (12) 



This latter result is exact at the bare one-loop level |T^ . 

The contribution of order g^ involves loop integrals with soft momenta, which requires 
using the HTL approximation for 11, where |16|,II3: 



fl^(^,A;) = ml,[l-tlog^], (13) 

UT{u;,k) = l[ml + {^-l)flL], (14) 

with rriD = gTJN/S. The spectral density of the corresponding gluon propagator consists of 
quasiparticle poles with momentum-dependent effective masses and Landau damping cuts 
for \uj\ < k. When k ^ gT, the additional pole associated to the collective longitudinal 



excitation has exponentially vanishing residue |17]. 

The order-gf'^ contribution in iShtl can be isolated as 

sSk/N, = /^i{2ReIl^[lm [u^ - k^ - fl^]- 

" 3 

J J 247r 

For the same reason as in the above scalar example, this is only part of the full order g^ 
contribution. Remarkably, it turns out that, as in the scalar case, Eq. ([15| ) is precisely 1/4 
of the correct result 5**^^^ = Ngm^^/(37f). 
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(a) (b) 

FIG. 2. NLO contributions to 6IIt at hard momentum. Thick dashed and wiggly lines with a 
blob represent HTL-resummed longitudinal and transverse propagators. 

The remaining order- (7^ correction comes from Eq. ([TT]) with ^IIt in place of Ilr, SUf^iy 
being evaluated at order gih'jj in HTL-resummed perturbation theory. The expression (|TTp 
is dominated by hard momenta k r^ T, and to the order of interest SUt is given by the two 
contributions shown in Fig. 2, in which one internal line is hard, and the other one is a soft 



resummed longitudinal (a) or transverse (b) gluon propagator. Diagram (a) restores the 
correct combinatorial factor of the longitudinal ring diagrams, whereas diagram (b) com- 
pensates for spurious transverse plasmon effects that are present in the HTL approximation 



SuTL 0. 



We turn now to the numerical evaluation of iShtl and shall discuss the effects of the above 
NLO contributions further below. iShtl involves two physically distinct contributions. One 
corresponds to the transverse and longitudinal gluonic quasiparticle poles, 



cQP _ _ /v 



Sk d 



Iik^wl'"°^^'-^-""'^ 



(27r)3 dT 

+^1^S l_e-./T I (16) 

where only the explicit T dependences are to be differentiated, and not those implicit in the 
HTL dispersion laws ujT{k) and ujL{k). The latter are given by the solutions of uj"^ — k"^ ~ 
tlTiyJTi k) and k"^ = —IIl{ul, k) 



Secondly, there are the Landau-damping contributions which read 

oo , 2 -.-. k 





-2 Im LIt Re [u;^ - A;^ - n^] "^ 
+ arg[fc2 + n^] - Imfli,Re[fc2 + n^]-i}. (17) 

The usual perturbative gf^- contribution (|ll]) is contained in the first term of Eq. (0); all 
the other terms in Eqs. ([T6|) , (p!^) are of order g^ in a small-(? expansion. 

In Fig. 3, we compare the numerical evaluation of iShtl/i^sb with the lattice data for 
purely gluonic QCD from Ref. [0, in the same manner as done in Ref. |T^, i.e., we use the 
two-loop running coupling constant ag^fi) of the MS-scheme with ^tto = Tc/1.03 and the 
renormalization scale is varied in the range /i = (^ . . . 2) x 27rT to give an estimate of the 
theoretical uncertainty. 

The thick dark-gray line represents the lattice data for the entropy density with the 
thickness of the line giving roughly the error reported in Ref. [|T^. Our result reproduces 
the lattice data rather well already for T > 2Tc. 
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1.5 2 2.5 3 3.5 4 4.5 t /T^ 
FIG. 3. HTL-improved results for the 2-loop entropy 5/5sb in purely gluonic QCD (full lines) 
with fJL varied between ttT and 47rT; our estimates for NLO effects are given by the dash-dotted 
lines. The lattice result for the entropy is represented by the dark-gray band. For comparison, the 
HTL-resummed results of Ref. |19| for the 1-loop pressure are given by the dotted lines, the lattice 
results for P/PsB by the light-gray band. [Because of the weak temperature-dependence of the 
theoretical results the predictions for 5/5sb are approximately those for P/PsB and vice versa.] 



Also given in the same figure are the lattice data for the pressure P/Psb (the lower light- 
gray band) and the result for the full HTL-resummation of the one-loop pressure reported in 
Ref. IQ (dotted lines). The discrepancy between the latter can be attributed in part to an 
overcounting in the interaction pressure in Ref. |l^, which will be corrected only in a fully 
resummed two-loop calculation. By contrast, our approach has the advantage of including 
the correct leading order interaction terms already in a pure HTL approximation, and also 
of manifest ultra-violet finiteness, thus avoiding the introduction of artificial counter-terms 
depending on the thermal mass. 

Turning now to the NLO approximation to S, we note that, at hard momenta k, 
Re6IlT{uj^ = k^) is, unlike m^, a nonlocal quantity. In order to get some estimate 
on its numerical effect, we approximate it by a constant correction which has the right 
magnitude to produce the known perturbative coefficient of order g^. We choose this as 
m^ + 5m^ = g'^T'^N /{6[1 + \/3Ng/n]), and include this correction only in the hard mo- 



mentum region defined hj k > M = y/2'n-Tmo- The reason for tliis is tliat we do not want 
to cliange by liand tlie overall scale of HTL contributions in the soft regime, where NLO 
contributions are known to behave quite differently: the long-wavelength plasma frequency 
receives much smaller negative corrections pOl, and the Debye screening mass is known to 



be even substantially increased [5^,53]. 



The dash-dotted lines in Fig. 3 give the correspondingly modified numerical results. In 
addition to a variation of /2 we have also included a variation of M, the boundary between 
hard and soft momenta, by a factor of 2 around its central value. These results happen to 
describe the lattice data surprisingly well, although their primary significance is to demon- 
strate the relative stability of our scheme upon inclusion of terms that restore equivalence 
with the known perturbative result up to and including order g^. A full NLO calculation 
still remains to be done and will be presented in a forthcoming publication. 

This work was supported by the Austrian-French scientific exchange program Amadeus 
of APAPE and OAD. 
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